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Abstract 

We study the energy level spacing of perturbed conformal minimal models in finite volume, 
considering perturbations of such models that are massive but not necessarily integrable. 
We compute their spectrum using a renormalization group improved truncated conformal 
spectrum approach. With this method we are able to study systems where more than 
40000 states are kept and where we determine the energies of the lowest several thousand 
eigenstates with high accuracy. We find, as expected, that the level spacing statistics of 
integrable perturbed minimal models are Poissonian while the statistics of non-integrable 
perturbations are GOE-like. However by varying the system size (and so controlling the 
positioning of the theory between its IR and UV limits) one can induce crossovers between 
the two statistical distributions. 



1 Introduction 



The distinction between integrable and non-integrable models is a central problem in many 
fields of physics and mathematics. The question is important in the quantum case, particu- 
larly for those systems with infinite degrees of freedom and local interactions. In this case, 
the only general result available consists in the identification of one-dimensional systems as 
the only models where integrability survives to the quantum level [T]. This no-go theorem, 
while limiting quantum integrability to one-dimensional theories, has, at the same time, led 
to intense focus on these same systems, both in their lattice and continuum formulations, 
leading to a number of remarkable developments in many-body physics (see, for instance 
[SJEIIHEIEIEIIEIIH]). The precise control offered by integrable field theories allows one 
to also examine non-integrable models, both in the massive [10] or massless cases jllj . 
through perturbation theory based on the exact matrix elements of the perturbed opera- 
tor. In such cases the nature of the non-integrable theory depends crucially on whether the 
perturbation is local or non-local with respect to the excitations of the original integrable 
model. In the former case, the mass of the excitations moves proportionally to the strength 
of the perturbation and therefore the excitations below the lowest threshold remain stable 
while those above threshold acquire an imaginary part in their mass and then decay \1'2\ . 
In the latter case, much more dramatic results occur. There the non-local nature of the 
perturbation gives rise to confinement phenomena of the excitations, with a consequent 
drastic change in the spectrum pm EJ H3] (see also [H H3 EH [UJ US Q2]). 

Despite an impressive mass of knowledge collected over the years following the Coleman- 
Mandula theorem, it is fair to say that we are still missing a well-established set of theo- 
retical tools for studying in a satisfactory way several aspects of non-integrable systems. 
Lacking, for instance, is a correspondent to the KAM theorem [20] that in classical me- 
chanics describes the effects of the breaking of an integrable Hamiltonian and gives the 
conditions under which quasi-periodic motion survives. Also lacking, significantly, is a 
proper understanding of time evolution of quantum systems with infinite degrees of free- 
dom. The importance of this issue and its intertwining with the system's integrability or 
non-integrability is seen in the series of experiments on one-dimensional cold atom sys- 
tems. Depending on the dynamics and the boundary conditions to which these have been 
subjected, the experiments have produced two sets of results. In one, the time evolution 
of the systems is insensitive to the initial conditions with a consequent thermalization of 
the final state |3T]. In the second |22j, the system has shown no sign of thermalization, i.e. 
the systems did not equilibrate even after thousands of collisions of the particles. 

The aim of this paper is to demonstrate a further discriminator between integrable and 
non-integrable quantum systems. Obviously, a general strategy in identifying integrable 
models is to show that the model admits a number of non-trivial conserved quantities 
equal to the number of degrees of freedom. This number becomes infinite in the case of 
a quantum field theory. Following this route is a difficult task for many models. There 
are, however, other features that provide hints of whether a given quantum system is 
integrable or not. The structure of energy level spacings in a quantum model is one such 
feature [27]. In particular the presence of an infinite number of conserved charges should 
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give rise to a high degeneracy in the energy levels, while a non-integrable system should 
present no sign of such degeneracies, provided we take into account "trivial" symmetries 
such as translational invariance, discrete symmetries, and so on. 

In the light of these remarks, we focus our attention here on the spectrum of a special 
class of quantum Hamiltonians obtained by deformations of the minimal unitary models 
of conformal field theories (CFT) [3J. The minimal unitary models, whose Hamiltonian 
is denoted by the sequence H p (p = 3,4, ■ ■ • ), are the continuum limit of certain (often 
elementary) quantum lattice models. The simplest of the sequence, H 3 , corresponds to 
the Ising model at its critical point. The next in the sequence, H4, is the tricritical- 
critical Ising model. is the three-state Potts model, and so on. Alternatively, the 
sequence corresponds to a set of Landau-Ginzburg theories where the potential of the order 
parameter field, $, is a polynomial of order $ 2 ( p_1 ) |25j . Or equally good, the quantum 
Hamiltonians H p are those which enter the transfer matrix of the two-dimensional RSOS 
lattice statistical models with maximal height p + 1 [26J. 

We focus on these theories for several reasons. First and foremost the minimal models 
and their perturbations offer a rich variety exhibiting all the different desired classes of 
behaviour. CFT's are themselves integrable field theories and for the minimal model sub- 
class, everything can be explicitly calculated. Furthermore perturbations of minimal CFTs 
can be both integrable jl] and non-integrable. Secondly, CFT's and their deformations, 
irrespective of whether they are integrable or not, admit a Hamiltonian formulation in a 
geometry of an infinite cylinder of circumference R. As shown by Al. Zamolodchikov |23| . 
the matrix elements of the Hamiltonian in this geometry can be obtained using conformal 
data alone. 

This formulation allows one to numerically study such theories through what is known 
as the truncated conformal space approach (TCSA). The "truncated" in the TCSA refers 
to the need to truncate the Hilbert space at some level (even in finite volume, the minimal 
models are defined on a continuum space and so their corresponding Hilbert space is infinite 
dimensional). While all aspects of the perturbed minimal models can be accessed in this 
fashion (despite the truncation) - including the masses of the lowest excitations, matrix 
elements, and even decay rates of unstable particles [24J - we focus here on the energy 
levels. Because we are interested in level spacing statistics (rather than, more typically, 
the particular value of one low-lying energy level), we need to accurately compute (at 
least) on the order of 4000-5000 levels of the system. In order to do this in the context 
of the TCSA, we employ a numerical renormalization group adapted to the particulars of 
the TCSA |45j . The numerical renormalization group used here is akin to the Wilsonian 
numerical renormalization group used in the study of quantum impurity problems. It 
is then a variational approach where the full Hamiltonian is approximated by a matrix 
product ansatz. To ensure that the best possible ansatz is employed, we introduce here 
a sweeping procedure. This ensures the thousands of needed energy levels are accurately 
determined. 

The paper is organized as follows. In Section 2, with the help of some simple examples, 
we remind the reader of the main motivations to study a model's level statistics. We also 
state the result of a general theorem that a set of uncorrelated variables has level spacing 
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statistics governed by a (generalized) Poisson distribution. In Section 3 we present both 
how the TCSA works and the necessary details of the conformal minimal models relevant 
to its implementation. In particular we describe how the numerical renormalization group 
is implemented with the TCSA and how we have expanded the procedure to include sweeps 
to optimize results. In Section 4 we present our results on the level spacing statistics for a 
number of integrable and non-integrable models. We also demonstrate that by varying the 
system size R (at fixed perturbation strength) we can induce a crossover from the Poisso- 
nian statistics applicable to the level spacings of an integrable model to the level spacings 
of a non-integrable model. This crossover occurs generically because at sufficiently small R, 
the perturbation is weak and the level spacing statistics is that of the unperturbed model, 
while at sufficiently large R the truncation effects of the TCSA render any perturbation 
non-integrable. 

2 Level spacing statistics in non-integrable and inte- 
grable models 

Michael Berry originated the idea of studying level spacing statistics to analyze integrable 
and non-integrable behaviour [27] . In particular he was interested in the question of how 
a quantum system (which is in a sense deterministic) can exhibit chaotic motion in the 
classical limit. In this respect, a quantum integrable system was conjectured to exhibit 
deterministic motion in this limit, while in this same limit a non-integrable quantum system 
sees chaotic motion. The spectrum of the families of quantum systems was then expected 
to show distinct features, in particular regarding the distribution of level spacings. 
Berry conjectured that the level spacing of an integrable system is Poissonian: 



while the spacing of a non-integrable system is a Wigner (or Gaussian Orthonormal En- 
semble, GOE) distribution: 



where s is the normalized spacing. If {E}f =1 is the set of energy levels of the systems, the 
normalized spacing corresponding to the i-th level is s, = — E^i) p(Ei) where p(Ei) 
is the density of states at energy Ei [28j |29| [37]. This normalization ensures the mean 
spacing is 1. 

While Berry himself provided good arguments for the integrable case [27J, to the best of 
our knowledge, there is no general proof of this conjecture although it has been supported by 
several subsequent studies (see, for instance, [271 129| [30j ETJ [321 E31 EH [35]). While dealing 
with any particular system, care needs to taken in calculating the spacing distribution. 
One must work in the subspaces where all the quantum numbers of the trivial symmetries 
that a system might have (e.g. translational invariance, discrete symmetries, spin reversal) 
have been fixed. In a recent paper [37] an unexpected spacing distribution was found for 
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the spin zero sector of a non-integrable XXZ chain. The spacing distribution resembled 
an average between a Poissonian and a GOE. Later the authors realized [38] that a spin 
reversal symmetry, present in the spin zero sector alone, was not taken into account. Taking 
into account such degeneracies resulted in a GOE distribution as expected. 



2.1 Some basic considerations and examples 

Before considering the level spacing of perturbed minimal models, we will examine a toy 
model so as to gain some intuition. Imagine a quantum Hamiltonian depending upon a 
parameter I. Plotting the energy levels vs. I, we should be able to see levels crossing one 
another if symmetries are present, while we should see instead avoided crossings in the 
absence of all symmetry. The reason behind this is that in a given symmetry subspace, 
states are coupled together by the Hamiltonian, while in different subspaces, states are 
decoupled and so evolve independently as we vary 




Figure 2.1: A set of lines with random slopes and intercepts. The spacing distribution is 
evaluated for a specific value of I, indicated by the red line in the figure. 



A naive way to model such a situation is to consider a set of straight lines on the 
plane, whose slope and intercept are independently identically distributed (i.i.d.) random 
numbers with a flat probability density (see Figure 2.1 ). Consider the y-axis as the energy 
and fix a certain value of I. The values of the ordinates of the lines at this I then are 
interpreted as a set of energy levels. The level spacing distribution of such a set of energy 
levels is shown in Figure 2.2. The distribution is clearly a Poissonian. This result could have 
been guessed from a remarkable theorem in probability theory by Pyke [36J. He showed 
that if we consider i.i.d. random numbers with probability density f(x), the limiting level 
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Figure 2.2: Spacing distribution for 100 000 random lines on a plane. The continuous line 
is the Poissonian distribution. 



spacing distribution is given by 



.'/(*) = / f( x ) e sf(x) dx; 
's 



f(x)dx = 1. 



(2.3) 
(2.4) 



where S is the support of f(x). For f(x) = const this results in a Poissonian distribution. 
It is interesting to notice that it is impossible to obtain a g(s) that tends to zero for s = 
(like the GOE) since this would imply J s f 2 (x)dx = 0, a condition clearly impossible. This 
is by no means unexpected. As we argued before, the levels in a non-integrable system in 
a given symmetry sector are correlated and cannot be mimicked by independent variables. 
On the other hand, the levels of an integrable system are split into an infinity of distinct 
independent subspaces and so behave largely as uncorrelated variables. 

We may think that the reason why a distribution with g(0) = is impossible to get 
from i.i.d. random variables is their independence. To test such a statement we generated 
a set of correlated Gaussian random variables, with mean value /i = 0, variance a 2 = 1 and 
correlation r. The trick to generate such a sample is straightforward. Start with N + 1 
independent Gaussian random variables Xi, .., x^, y with fi = and a 2 = 1. Define the 
linear combinations 



XiVT 

N 



C 



E 



r + yy/r 



(2.5) 



y/l + r(N-l) y/l + r(N-lY 
with < r < 1. Since these are linear combinations, the Zj, ( are still Gaussian random 



variables. However they are now correlated. Indeed, from definitions (2.5) 

(*) = 0; 
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(ZiZj) 
(ZiZj) 

(ZiC) 

(C 2 ) 



r, 
1- 
0; 



(2.6) 



as required. The variable ( is independent. Focusing on Zi, we obtain N Gaussian corre- 
lated variables with the required properties. The spacing distribution for a sample of 100 



000 of such variables is presented in Figure 2J3 for several values of r. Independent of this 
r, the distribution is still Poissonian. 

One may question whether a linear transformation is sufficiently complex enough to 
alter the spacing distribution. We thus try instead a quadratic transformation 



X; 



Vgj -r) + 2y/4^T 
6 ' 



(2.7) 



with i < r < 1. The new variables are no longer Gaussian, but have correlation r and 



variance 1. But again the level spacing distribution is Poissonian (see Figure 2.4) 



We conclude from this exercise that the correlation between variables required to get 
a GOE-like curve is much more involved. In random matrix theory the starting point is a 
random Hamiltonian H, symmetric for the GOE case, and so having JV ( j ^ +1 ) independent 
random variables. The joint probability distribution of the N eigenvalues corresponding 
to this H is given by 



P(E 1 ,....E N )ocl[(E j 

j<k 



E k )e 



1 p2 

2 2^i=l 



(2.8) 



The term Ylj<k(^j ~ * s ^ ne responsible for the peculiar level spacing correlation. It 
comes from the Jacobian of the transformation from the original matrix elements Hy to 
the energies Ei (for details see [46j). The two-point correlation function is then a function 
of the distance Ei — Ej and goes to zero for Ei 



Ej. 



Such a correlation is due to the 





(b) 



(c) 



Figure 2.3: Level spacing distribution for Gaussian correlated random variables with 100 
000 samples taken, (a) r=0.1; (b) r=0.5; (c) r=0.9. For comparison we present the Pois- 
sonian distribution as a continuous line while the GOE is given by a dashed line. 
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Figure 2.4: Level spacing distribution for non-Gaussian correlated random variables with 
100 000 samples taken and r=0.5. Continuous line: Poissonian distribution; dashed line: 
GOE. 

highly nontrivial operation of diagonalizing H, a procedure not mimicked properly in the 
examples above. 

We thus conclude that the energy level spacings of an integrable model behave as 
independent variables inasmuch as they satisfy Poissonian statistics. However as we have 
attempted to demonstrate, the energy levels of a non-integrable model must exhibit highly 
non-trivial correlations in order to give GOE statistics. 



3 Studying perturbed minimal models using the trun- 
cated conformal space approach 

In this section we discuss how the level spacing statistics of the perturbed minimal models 
can be studied using the truncated conformal spectrum approach. We begin by reviewing 
briefly some necessary details of the minimal models of conformal field theory. 

3.1 The conformal minimal models 

The minimal models are a sequence of two dimensional conformally invariant quantum 
field theories. They are labelled, M M , where p and q are two coprime integers (a unitary 
model is obtained for q=p+l). They have central charge of less than 1 given by 

c=l-6 {p - q)2 . (3.1) 
pq 

The field content of M VA consists of a set of primary fields, r>s , r>s with chiral dimensions 

K,s = h r , 3 = ^ - qS \ - (P - ?) , l<r<g, l<s<p. (3.2) 

4pq 
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The Hilbert space of these theories consists of highest weight states 



\h r , s ,h r , s ) = riS (O)0 rjS (O)|O), 

(where we take h = h so restricting ourselves to the diagonal theories), acted upon by 
combinations of Virasoro generators, L n ,L n , satisfying the algebra 

[L n , L m ) = (n - m)L n+m + 5 n+mfi —n(n 2 - 1). (3.3) 

The set of states obtained by acting upon a single highest weight state with the Virasoro 
generators is known as a Verma module. On a cylinder of circumference, R, where time is 
taken along the cylinder's axis, the Hamiltonian and momentum are given by 



2ir - c 2ir 

H p , q = ~^(L + L - —) , P m = — (L - L ). (3.4) 



(we denote the sequence of unitary models by H p and P p ). This set of Hamiltonians has 
an infinite set of conserved charges which can be written as normal ordered products of 
polynomials of the stress energy tensor |41] . 

We study deformations of the minimal models arrived at by perturbing them by some 
combination of fields in the theory: 

f R 

i Jo 

These perturbations in general destroy the infinite set of conserved quantities present at the 
conformal point. The only deformations which generically give rise to integrable theories 
away from criticality are perturbations involving one (and only one) of the fields 0i2, 
013 and 02i [!]• Perturbations involving multiple fields generically lead to non-integrable 
Hamiltonians [8| 1 15 j . 

3.2 Minimal form of TCSA 



To study the Hamiltonians in Eqn.(3.5) we employ the truncated spectrum approach. 
Because we have complete control over the conformal Hilbert space, we can enumerate 
all the states in the Hilbert space. Although the conformal level spacing is finite (and 
essentially equal to 27r/R), the continuum nature of the theory means the Hilbert space 
is nonetheless infinite dimensional. In order to make the theory tractable we truncate 
the spectrum, keeping states whose energy in the conformal field theory is less than some 
cutoff, E c . Being left with a finite Hilbert space, we can form the Hamiltonian matrix 
through the evaluation of the following set of matrix elements: 

(^\H\V) = ^(A* + A^-^)^ + ^^(*|0 ri!S .(O)|* / )^,s^, (3.6) 



9 



where |\&) and are two typical states in the truncated spectrum. Here Lq\^) = A^\ty) 
and we suppose the state has spin s. The matrix elements (^\(j) r . tS .(0)\^f') are given 

in terms of the structure constants C^ r ^' of the theory and the chiral dimension of the 
perturbing operator: 

It i' i 

where M^^i is some number depending on the particular Virasoro modes acting on the 
highest weight state in creating the state \^f), \ ^'). The structure constants of minimal 
models were computed by Dotsenko and Fateev [44J. Once these matrix elements are 
evaluated, we can easily diagonalize numerically the corresponding matrix and thus extract 
both the spectrum and matrix elements of the full (truncated) theory. 

For large values of the truncation energy, E c , computing this matrix can involve the 
computation of millions of matrix elements. As such, this computation has been automated. 
To reduce the numerical burden further we perform the computations in the independent 
subspaces of the Hamiltonian (typically a given perturbation of a minimal model will 
not mix together all the Verma modules). We also perform the determination of the 
independent subspaces automatically through the computation of the structure constants, 
Cy^ 1 . All told, the code we have developed to study the level spacings in perturbed 
minimal models can do so for arbitrary M Pi? and for arbitrary perturbations. 



3.3 NRG improvements on TCSA 

In performing this computation upon the truncated Hamiltonian, the larger we choose the 
cutoff, E c , the better. However increasing E c comes with a numerical cost. The number 
of states in the truncated theory grows exponentially with E c . For a given Verma module, 
the number of states will grow as roughly the number, P(N), of partitions of the integer 
N = RE c /2n: 

P(N) ~ — *-= e^v 7 ^, iV ^ oo. (3.8) 

For the minimal CFTs, this growth is slightly slower due to the presence of zero norm 
states, but nonetheless is still exponential. 

A way in which to improve upon the precision with which the high energy levels are 
handled without exponentially increasing the computational burden is given by a numerical 
renormalization group (NRG) approach. This approach, modelled after Wilson's numerical 
renormalization group for quantum impurities, allows one to keep hundreds of thousands 
of states in the truncated Hilbert states (as opposed to order ten thousand without NRG 

The basic idea behind the approach is the ability to order the states of the CFT in 
terms of their importance for determining the (relatively) low energy behaviour of the full 
theory. This ordering is determined by the unperturbed energies of the states of the CFT 
field theory. Because the perturbation of the CFT that we are considering is relevant, low 



10 



energy states of the CFT are more important for determining the low energy physics than 
high energy states are. The situation is analogous to Wilson's treatment of the Kondo 
problem. There Wilson mapped a Kondo impurity interacting with free electrons onto 
a tight binding model of non-interacting electrons on a lattice extending from to oo 
interacting with an impurity located at x — 0. Wilson reasoned that a site on the lattice 
close to the impurity will influence the physics more than a site further away. Thus the 
natural way to proceed in such a problem is to first diagonalize the interactions of the 
Kondo impurity with the electrons on nearby lattice sites, and then take into account in 
a progressive fashion, sites further and further away from the impurity. We mimic this 
procedure here. However instead of distance from the Kondo impurity as a metric of 
importance, we use the magnitude of the unperturbed energy of a conformal state. 

The NRG algorithm by which a truncated perturbed conformal field theory is treated 
is as follows. We imagine ordering in energy the states of the unperturbed conformal field 
theory, |0), |1), • • • . As a first step we take the N + A of the states lowest in energy. We 
form the full Hamiltonian (both CFT and perturbation) in this truncated basis of iV + A 
states. Just as with the ordinary TCSA, we diagonalize (numerically) the problem leaving 
us with N + A eigenstates, E\, ■ ■ ■ , E 1 N+A and eigenvectors, \E)\, • • • , \E)]^ +A . In the next 
step of the NRG algorithm, we order these N + A states in ascending order of their energies 
and toss away the top A states. We express these remaining eigenvectors in terms of the 
conformal basis as follows: 

\E)\ = A\ m \m) + B l kml \m' + N); k = l,---,N, (3.9) 

where A 1 is an N x N matrix and B l is an N x A matrix and repeated indices are 
summed on. To these N eigenstates we add the next A states from the unperturbed 
theory,|m + A + 1), • • • , \m + 2A). This leaves us again with a truncated Hilbert space 
of iV + A states. We reform the Hamiltonian in this new basis and then rediagonalize so 
extracting a new set of N + A energies, Ef, • • • , Efj +A , and eigenstates, \E)\, • • • , \E) 2 N+A . 
We again order the eigenstates in energy and toss away the top most A states. If we 
reexpress the remaining eigenstates, \E)\ in terms of the original conformal basis we obtain 

\E)\ = AlJE^ + BlJm + N + A) 

= (A 2 A 1 ) km \m) + (A 2 B 1 ) km/ \m' + N)+B 2 mt \m' + N + A), (3.10) 

where again A 2 is an N x N matrix and B 2 is an iV x A matrix. Here the index m runs 
from 1 to N while the index m' runs from 1 to A. 

We can repeat this procedure ad libitum: we order the new set of N + A states, toss 
away the top most A states, add the next A set of unperturbed states from the CFT, 
reform the Hamiltonian, rediagonalize, etc. In this way, we allow the higher energy states 
of the unperturbed CFT to mix in with eigenstates of the full (but truncated) theory. As 
we keep iV + A fixed at each step, the associated computational problem grows only as the 
square of the total number of states kept (owing to the need to manipulate states (albeit 
only iV + A of them) which are expressed in terms of an ever growing basis as the NRG 
proceeds). 
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At the nth-iteration, the eigenstates have the form 

\E) n k = (A n ---A 1 ) km \m) + (A n ---A 2 B 1 ) km r\m + N) 
+(A n --- A 3 B 2 ) km , \rri + N + A) + • • • 

+ (A n B n - 1 ) km ,\m' + N+(n- 2)A) + B n krn ,\m' + N + (n - 1)A). (3.11) 

We see that each term in the above sum has a matrix product state form. The approx- 
imation encoded in the NRG is then one where we study a Hamiltonian arrived at by 
projecting the original Hamiltonian onto a space composed of matrix product states of the 
above form. 

The NRG scheme as described above is able to obtain accurately the low energy portion 
of the spectrum. However as we interested in determine the level spacing distribution we 
need the first several thousand levels to high accuracy. In order to obtain such accuracy 
we adopt a sweeping procedure not dissimilar to the finite volume algorithm in DMRG 
|39j. Or equally good, this procedure may be considered as an extended set of Jacobi 
transformations done to diagonalize a symmetric matrix (extended in the sense that we 
are zeroing blocks not individual elements of a matrix). The procedure works as follows. 
Suppose we have completed the NRG procedure as described above for M iterations. And 
in doing so we would have run through the first N + MA states of the conformal basis. We 
now begin anew. But instead of working with N + MA conformal states, we work with a 
basis formed from the eigenstates generated in the the NRG procedure. This basis is given 
by 

\ E )'i > ■■■■> \ e )ni \ e )n+i, ■ ■ ■ > \ e )n+ai \ e )n+h ■ ■ • > \E)n+a> ■ ■ ■ 

\E) N +i, • • • , \ e )n+L \ e )n+h • • • j \ e )n+a- (3.12) 

We see that this basis is formed from taking the first N states coming from the last NRG 
iteration (and so the best guess we have at the N lowest energy states in the theory), 
followed by the M x A states we discarded in the M iterations of the NRG. This basis 
is much closer to the true eigenstates of the system than the initial conformal basis was. 
Using this basis we repeat the above described set of NRG iterations (although because 
this basis is not orthonormal an additional step at each iteration is needed to form an 
orthonormal basis). This set of iterations constitutes a single sweep. We find that the 
eigenenergies converge rapidly after only a handful of sweeps. To test this procedure we 
compared exact diagonalization of large matrices (12000 x 12000) with the NRG plus three 
sweeps with iV = 4000 and A = 1000. We found that this number of sweeps was enough 
to guarantee that 90% of the first N eigenstates matched the exact diagonalization to 4 
significant digits. 
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4 Level Spacing Results 



The results of any study of the level spacings of a perturbed minimal model will depend 
crucially on the value of the system size, R. For a generic perturbation, the model will be 
massive with some mass scale m = This allows us to identify two different regimes: 
R <C £ and R ^> £. In the regime R <C £, the theory is in its UV limit. Regardless of 
whether the perturbation is integrable or otherwise, the level spacing distribution will be 
close to being Poissonian. On the other hand for R ^> £, the theory is in its IR limit 
and whether the theory is integrable or non-integrable becomes important for determining 
the distribution its level spacing obeys. A consequence of this is that as R is varied for a 
non-integrable model the level spacing distribution will evolve from Poissonian to GOE. 




12 3 4 5 1 2 3 4 5 

(a) " (b) 

Figure 4.1: Energy levels(multiplied by R) vs cylinder size, (a) Integrable case, showing 
several level crossings, (b) Non-integrable case, showing instead avoided crossings. 

The truncation energy, E c used in studying a given theory with the TCSA also influences 
the level spacing distribution. To see this first imagine that we are in the IR limit and we 
plot the energy levels as a function of R for the integrable and non-integrable case. Further 
imagine that the truncation is taken to infinity (there is no truncation of the spectrum). 
In the integrable case, we see that energy levels cross at various values of R (schematically 
pictured in Figure 4.1a). This occurs because even though we work in a given symmetry 
sector where all the quantum numbers of the trivial symmetries are fixed, the eigenstates 
still carry the quantum numbers of the remaining non-trivial symmetries that are present 
in the integrable model. If these quantum numbers are different, the eigenstates cannot 
mix and their corresponding energies may cross. This crossing leads to weight in the 
level spacing distribution for arbitrarily small level spacings. While it does not explain 
the appearance of a Poissonian distribution, it is at least consistent with it. If however 
integrability is not present, eigenstates will be able to mix generically and what we will see 
instead is a set of avoided crossings (see Figure 4.1b). 

Now imagine that we include the effects of E c . The presence of a finite truncation of 
an integrable theory will generically destroy integrability. (The infinite set of conserved 
charges in the theory, being polynomials of the stress-energy tensor, involve Virasoro modes 
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of arbitrary high grade. But in the presence of a finite E c , modes with grade greater than 
RE c /27r are discarded from the theory so breaking the charges' structure.) Correspondingly 
level crossings will turn into avoided crossings. Thus a theory with a finite E c should not 
have a perfect Poissonian distribution. But of course if E c is chosen to be large enough, the 
breaking of integrability will be sufficiently weak that the distribution will remain nearly 
Poissonian. One reason it is necessary to use the numerical renormalization group in this 
case is to be able to set E c high enough so that this breaking is in fact small. 

The effects of E c are not independent of the system size, R, the relevant dimensionless 
parameter being ^E c /2tt = C,2n c /R, where n c is the level at which we truncate the chiral 
Verma modules. Thus for even a large value of E c , a sufficiently large value of R will break 
integrability. Thus for an integrable model, as we vary R we expect to see the level spacing 
distribution to evolve from Poissonian to something GOE-like. In the integrable case, we 
will be particularly interested in values of R large enough such the theory is in its massive 
phase (i.e. R > £) but where £E c /2ir is sufficiently large that integrability breaking is 
weak. 

To quantitatively characterize a given distribution g(s), it is common (see [50J for 
example) to define a quantity 77: 

= f™ ds\g(s) - e-\ 

11 JJ» da |H e -f- a _ e -.|- l4 - lj 

Note that our definition of r\ is differs slightly from the common one. Usually the integration 
is restricted in the region < s < Sq, where sq is the first intersection point between 
Poissonian and GOE. Such a definition completely disregards the tail of the distribution. 
However the tail is less sensitive to the bin size chosen, and has very different behavior in 
the Poissonian and GOE case. Therefore, disregarding the tail gives a much less reliable 
estimate for 77. 

From the above definition, we can extract the limiting behavior of rj. If rj is close to 1, 
g(s) is GOE-like. If instead it is close to 0, g(s) is Poissonian. And finally if it is close to 
|, the distribution g(s) is similar to (Poisson+GOE)/2. 



4.1 Tricritical Ising model 

We now consider as a first example the level spacing statistics of the perturbed tricritical 
Ising modeQ (for previous TCSA studies of this model see [T4"t EE]). The tricritical Ising 
model (TIM) is described by the minimal model M 4)5 . The field content is characterized 
by six primary fields. We will consider the perturbations by the fields 0^2 (the leading 
energy perturbation) and 02,2 (the leading spin perturbation) and a combination thereof. 
The scaling dimensions of these two fields are respectively ^ and Jjj. The perturbation 
02.2 is odd under Z2 symmetry while 01,2 is even. In analyzing these theories, we focus on 

■""We have also performed similar calculations, with similar results, for the usual Ising model, but we 
prefer to focus our discussion here on the TIM for its larger number of integrable and non-integrable 
deformations. 
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the zero momentum subsector. In the cases where o — > —o is a good symmetry [a is the 
spin field) we analyze the largest subspace. 

The first results we present are level spacings obtained for the integrable perturbation 
<pi 2 of the TIM. The coupling constant is set to 0.092834 . . ., to ensure a mass gap, m = 1 
[47]. For a truncation at level n c = 14, the zero momentum Hilbert space splits in two 
subspaces of 22559 and 18751 states. We employ the NRG with a base matrix of iV = 
4000, a step size of A = 1000, and use three sweeps to obtain convergence. Of the 4000 
eigenvalues in the base matrix, we employ the lowest 3500 in computing the level spacing 
distribution. For smaller values of R we are in the conformal limit. 




Figure 4.2: Level spacing distribution for M 4 5 + 12 for three values of mR: (a) mR=18.9, 
77=0.25, correlation coefficient for the fit to a Poissonian distribution, rp = 0.98; (b) 
mR=22, 77=0.22, rp = 0.99 (c) mR=25.2, 77=0.55, rp = 0.95, correlation coefficient for 
the fit to (Poissonian+GOE)/2, rAv — 0.98. In these three plots we also display a pure 
Poissonian distribution (solid line) and a pure GOE distribution (dashed line). 

In the pure conformal limit, all energy levels are of the form ^(n — 1/12) where n is 
an integer. Moreover, the degeneracy of each level n grows exponentially in n. Thus we 
might expect a disproportionate weight at zero in the energy level spacings. This is what 
we find in Figure 4.2a. 

As we increase the value of R we see that we obtain a level distribution that is Poissonian 
at all values of s (see Figure 4.2b). As we increase R further we see that the Poissonian 
distribution evolves into a distribution closer to a mix of Poissonian and GOE (Figure 
4.2c). As we have argued, this is to be expected as for larger values of R, we expect the 
integrability of the theory to be compromised by the presence of a finite energy cutoff. 

For the next example, we consider the non-integrable 2 ,2 perturbation of the TIM. 
Again we chose the coupling constant so that m = 1. The truncated Hilbert space consists 
of 41310 states. The NRG base matrix N, step A and number of states considered are 
chosen as in the <f)i^ case. The level spacing distributions for various values of mR are 
shown in Figure 4.3. For very small values of mR we see that the distribution is Poissonian 
as expected. As mR is increased the distribution becomes more and more GOE-like. We 
may ask whether this effect is aided by truncation effects. We see to some degree that it is. 
At the values of mR shown, the distribution is most GOE like for the largest value of mR 
(25.2). But from our study of M^^+tpi^ we know that at this value of mR the Poissonian 
distribution due to the model's integrability has already begun to break down. 
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Figure 4.3: Level spacing distributions for M 4j5 + 02,2 (a) mR=8.2, rj = 0.25, rp = 0.98; (b) 
mR=13.2, rj = 0.52, tgoe = 0.91 (the correlation coefficient w.r.t. the GOE distribution), 
r AV = 0.99; (c) mR=18.9, r) = 0.64, r GOE = 0.95, r AV = 0.98; (d) mR=22, r? = 0.68, 
tgoe = 0.97, r AV = 0.98; (e) mR=25.2, rj = 0.75, r GOE = 0.98. 



As a final case for the TIM, we consider perturbing the theory by both 1)2 and 02,2 
simultaneously. The coupling constants are chosen such that m — 1. The Hilbert space 
and N and A have the same size as the 02,2 case. The results are presented in Figure 4.4. 
As with + 02,2, we see that for + 0^2 + 02.2 increasing the value of R leads to a 
distribution that is more and more GOE like. 



4.2 Tetracritical Ising Model 

To demonstrate that our level spacing results possess a certain universality, we consider 
the tetracritical Ising model, a theory described by the minimal model M 5 6 . (Its partition 
function differs from three state Potts in that it is diagonal in the Virasoro characters.) 
The field content of this theory is characterized by ten primary fields. We will consider, 
separately, the perturbations 0^2 (integrable) and 02,2 (non-integrable). The scaling di- 
mensions of these fields are respectively | and ^. 

We take the truncation level to be n c = 11 with the corresponding size of the truncated 
Hilbert space being 27931 states for 2 2 case and 17601 for 12 . We set the coupling 



constants in both cases such that the mass gap is again 1. We see in Figure 4.5 that, 
as expected, the integrable perturbation has Poissonian level spacing statistics while the 
non-integrable perturbation has GOE statistics. 
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Figure 4.4: Level spacing distribution for + + 02,2 (a) mR=13.2, 77 = 0.60, r^os = 
0.94, r Ay = 0.99; (b) mR=18.9, r] = 0.74, r GOE = 0-98; (c) mR=22, 77 = 0.79, r GO £ = 0.99; 
(d) mR=25.2 77 = 0.80, r GO £ = 0.99. 

5 Conclusions 

In this work we have studied the level spacing statistics of perturbations of the conformal 
minimal models. We find that the level spacing statistics of the integrable variants of such 
models are Poissonian while their non-integrable counterparts obey GOE statistics. This 
general result does come with two caveats. For values of the system size, R, of up to 10 
times the correlation length, £, the system exhibits Poissonian-like statistics regardless of 
the whether the perturbation is integrable or not, a result of being in the conformal limit 
(at least at higher energies) for moderate values of R. The second caveat is that at very 
large values of R, values where ^E c is less than 7, the level spacing distribution becomes 
GOE-like regardless of the presence of integrability. This results from the presence of a 
finite energy cutoff in our numerical treatment of these theories. At large values of R (or 
equivalently small values of the dimensionless parameter E c t; = ), the presence of this 
truncation leads to a strong breaking of any symmetries in the model (and so integrability), 
so guaranteeing GOE like statistics. 

To determine the level spacing statistics we employed the truncated conformal spectrum 
approach. In this method the underlying conformal spectrum is truncated at some energy 
E c . In order to obtain results that approximate the theory without cutoff, we needed to 
chose E c to be large enough that the considered Hilbert space contained in certain cases on 
the order of 40000 states. This cannot be handled directly by exact diagonalization. Thus 
we used a numerical renormalization group procedure jl5] which enables large truncated 
Hilbert spaces to be studied in a numerically feasible fashion. While this approach yields 
accurate results for low lying eigenstates, it, by itself, cannot guarantee accuracy over the 
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Figure 4.5: (a) Level spacing distribution for M^ t % + 01^; mR=18.9, 77 = 0.20, rp = 0.98; 
(b) Level spacing distribution for M 5i6 + 0i i2 ; mR=18.9, r\ = 0.75, tgoe = 0.98. 

thousands of levels needed to determine the level spacing distribution. Thus we developed 
a sweeping procedure guaranteeing convergence of high energy eigenstates to their true 
value. 

In this paper we report the results of our analysis for two particular minimal models, the 
tricritical Ising model and the tetracritical Ising model. These results point to the infinite 
set of quantum Hamiltonians associated to the deformations of conformal minimal models 
as an ideal playground for studying important characteristics of generic non-integrable 
quantum systems. It would be interesting, for instance, to employ these models in order 
to test the validity of the various scenarios put forward for the thermalization of generic 
isolated quantum systems following a sudden quench [49J. Moreover, properly modified, 
the numerical method developed here can be also used to address the computation of 
non-equilibrium and time-dependent correlation functions of various order parameters. 
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